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FORMULATION OF ANALYSIS PROBLEM 


For an axisymmetric body Immersed in a propfan slipstream, or a jet, the 
effect of freestream nonuniformity must be accounted for to calculate the wave 
drag In a transonic flow. Since the flow field is rotational, in the present 
approach a rotation function (F) and a velocity function (<J>) are introduced in 
Euler's equations to result in a governing equation which is similar, in 
mathematical structure, to a full-potential equation as shown in figure 1 
(refs. 1 and 2). The equation is solved with the algorithm of reference 3. 
Following reference 3, the equation is cast in curvilinear coordinate systems 
with a body-normal coordinate system covering the front portion and a sheared 
cylindrical system used in the aft. The rotation function is calculated 
through Crocco's relation. In figure 1, only the equation in body-normal 
coordinates is shown. Details can be found in reference 2. 


Equation in Body Normal Coordinates: 
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• Velocity Components: 

u » g <t>^ + ( 1 + F ) cosG 


v = <t> n — ( 1 +F)sin0 


• Crocco's Relation for the Rotation Function: 
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T o = Stagnation temperature 


P 0 = Stagnation pressure 
Figure 1 



1086 



DESIGN FORMULATION 


The analysis program is coupled with the CONMIN optimizer (ref. 4) to 
design an axisymmetric body for minimal wave drag. The constraints include a 
specified maximum thickness and a tail thickness. To reduce the number of 
design variables, the shape is represented by a modified Fourier series with 
its coefficients being the design variables (fig. 2). 


• Objective Function to be Minimized: 

OBJ - -0.1/(0.001 + C u ) 

• Constraints: 

G(2)-10(1-r m „/r t )<0 

Q < 3 >- <° 

<3(4) - 1 - r J r, <0 

where r u , ij, t u , and t 4 are the specified upper and lower bounds of maximum thickness 
and the tail thickness, respectively 


Body Shape Representation: 

1 . Rounded Nose and Tail 


r = |/^i(9 + sine)-%i 


(6 - Sine) + ^L(e - + £ A n . single j! 


2. Rounded Nose only 


r = 2 |an+i sin0 cos|- + ^ A n cos(n-1)e|, 9 = cos’ 1 |y- -1 j , i= body length 


Figure 2 
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DESIGN ALGORITHM 


The design process is started by identifying the design variables from 
the input shape through a Fourier analysis. By perturbing the design 
variables (i.e. , the Fourier coefficients) one at a time, gradients of the 
objective function and constraint equations can be calculated. To reduce the 
computing time, these gradient calculations are made with a small change in 
the design variables. Typically, the change (to be called the step size) is 
taken to be 0.1% ~ 0.5% of each design variable, but not less than 
0.00035 ~ 0.0005. These gradients are all calculated with the same starting 
Rvalues. These <|>-values are updated if the design is feasible (fig. 3). 



Figure 3 
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SUMMARY OF SOME OPERATIONAL EXPERIENCE 


1. To reduce the computing time, an improved design was chosen by comparing 
all calculated results at the same numerical accuracy. That is, in 
analysis, convergence to a very small residual was not demanded. 

Instead, evaluation of all gradients was based on convergence to the same 
maximum residual. All starting values of the velocity function (•)>) for 
gradient calculation were the same. The velocity function was updated 
only when the design has improved. 

2. Typically, 2 to 3 iterations were performed in each run to allow manual 
adjustment in step size for gradient evaluation. If the design was not 
improved, the step size should be reduced. For this purpose, the 
solution was always saved in a file for possible re-use. 

3. Since a smooth input shape was highly desirable, it was found 
advantageous to Fourier-analyze the input shape separately and then use 
the resulting Fourier coefficients, or modified values if desired, to 
generate a starting shape with more defining coordinate points. 

4. All design exercises have been achieved with 81 x 81 grid points. 

Attempt with 41 x 41 grid points has not been successful. 
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DESIGN OF AN AXISYMMETRIC BODY WITH A FINENESS RATIO OF 8.33 

IN A NONUNIFORM FLOW 

The initial shape was assumed to be given by the NACA-0012 contour. Six 
design variables (A n ) were used. The tail thickness was constrained to be 
between 0 and 1%. The Mach number in the external flow was 0.98, and that 
over the body was 0.995. The results in figure 4 indicated that reducing the 
nose radius and increasing the thickness in the aft portion would reduce the 

wave drag (cj ) by 29%. 
a w 



Figure 4 
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DESIGN OF AN AXISYMMETRIC BODY WITH A FINENESS RATIO OF 5.0 
IN A TRANSONIC UNIFORM FLOW 

The initial shape was generated from the NACA-0020 contour. The 
freestream Mach number was 0.925. The step size used in gradient evaluation 
was 0.1% of the design variables with a minimum change of 0.0005. The results 
in figure 5 showed that to reduce the wave drag, the pressure peak in the nose 
region must be reduced. As a result, the shock strength could also be 
decreased. A reduction in wave drag by 41% was achieved. 



x/E 

Figure 5 


1091 


DESIGN OF AN AXISYMMETRIC BODY WITH A FINENESS RATIO OF 5.0 
IN A TRANSONIC UNIFORM FLOW 

The final body shape given in figure 6 indicated, again, that to reduce 
the wave drag, the nose radius should be reduced and the maximum thickness 
location moved aft. Further improvement could be made only if a larger tail 
thickness was allowed. 
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DESIGN OF AN AXISYMMETRIC BODY WITH A FINENESS RATIO OF 5.0 
IN A TRANSONIC NONUNIFORM FLOW 


The initial shape was again generated with the NACA-0020 contour. The 
freestream Mach number varied from 0.90 away from the body and a 0.95 near the 
body, with an average of 0.925. The step sizes in gradient evaluation used 
ranged from 0.1% at the beginning to 1%. The final value used was 0.5%. It 
was found that if the step size was greater than 0.5%, little improvement 
could be made. The results in Figure 7 showed that the change in shape 
successfully reduced the shock strength. The wave drag was reduced by 65%. A 
larger drag reduction was possible in this case perhaps because the external 
Mach number was lower than that in figure 5. 
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DESIGN OF AN AXISYMMETRIC BODY WITH A FINENESS RATIO OF 5.0 
IN A TRANSONIC NONUNIFORM FLOW 

Again, a favorable shape was one with a reduced nose radius and a thicker 
aft portion as shown in figure 8. Further change was difficult because of the 
constraints of maximum thickness and tail radius. 



1094 



CONCLUSIONS 


An inviscld transonic code capable of designing an axisymmetric body in a 
uniform or nonuniform flow was developed. The design was achieved by direct 
optimization by coupling an analysis code with an optimizer. 

Design examples were provided for axisymmetric bodies with fineness 
ratios of 8.33 and 5 at different Mach numbers. It was shown that by reducing 
the nose radius and Increasing the afterbody thickness of initial shapes 
obtained from symmetric NACA four-digit airfoil contours, wave drag could be 
reduced by 29% for a body of fineness ratio 8.33 in a nonuniform transonic 
flow of M ■ 0.98 to 0.995. The reduction was 41% for a body of fineness ratio 
5 in a uniform transonic flow of M - 0.925 and 65% for the same body but in a 
nonuniforra transonic flow of M ■ 0.90 to 0.95. 
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